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ABSTRACT 



Aims. We aim to demonstrate the usability of mm-wavelength imaging data obtained from the APEX-SZ bolometer array to derive 
the radial temperature profile of the hot intra-cluster gas out to radius 7-500 and beyond. The goal is to study the physical properties of 
the intra-cluster gas by using a non-parametric de-projection method that is, aside from the assumption of spherical symmetry, free 
from modeling bias. 

Methods. We use publicly available X-ray spectroscopic-imaging data in the 0.7-2 keV energy band from the XMM-Newton obser- 
vatory and our Sunyaev-Zel'dovich Effect (SZE) imaging data from the APEX-SZ experiment at 150 GHz to de-project the density 
and temperature profiles for a well-studied relaxed cluster, Abell 2204. We derive the gas density, temperature and entropy profiles 
assuming spherical symmetry, and obtain the total mass profile under the assumption of hydrostatic equilibrium. For comparison 
with X-ray spectroscopic temperature models, a re-analysis of recent Chandra observation is done with the latest calibration updates. 
We compare the results with that from an unrelaxed cluster, Abell 2163, to illustrate some differences between relaxed and merging 
systems. 

Results. Using the non-parametric modeling, we demonstrate a decrease of gas temperature in the cluster outskirts, and also measure 
gas entropy profiles, both of which are done for the first time independently of X-ray spectroscopy using the SZE and X-ray imaging 
data. The gas entropy measurement in the central 100 kpc shows the usability of APEX-SZ data for inferring cluster dynamical states 
with this method. The contribution of the SZE systematic uncertainties in measuring T e at large radii is shown to be small compared 
to XMM-Newton and Chandra systematic spectroscopic errors. The total mass profile obtained using the hydrostatic equilibrium as- 
sumption is in agreement with the published X-ray and weak lensing results; the upper limit on M 2 oo derived from the non-parametric 
' method is consistent with the NFW model prediction from weak lensing analysis. 

Key words. Galaxies: clusters: individual: Abell 2204 - Cosmology: observations - intergalactic medium - cosmic microwave 
background - X-rays: galaxies: clusters 
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■ 1. Introduction Accurately determining the thermodynamic state of the 

intra-cluster medium (ICM) out to a large radius is critical 

Current cosmological models are built upon two complementary for underst anding the link between cluster mass and observ- 

approaches of astronomical observation: the statistical study of aWes For oyer a decad6i observations of the thermal Sunyaev- 

the ensemble properties in a large sample of objects (i.e. from Zerdovich Effect (tSZE> here after simply SZE; Sunyaev & 

surveys) and the detailed analysis of the individual objects for Zel > dovich 1970> Birkinshaw 1999) have been considered as a 

gaining better understanding of the physical processes affecting promising comp i e ment to X-ray observations for modeling the 

those ensemble properties. This is particularly important in the ICM in galaxy Qnly recently hafj k been pQSsible 

study of galaxy clusters, where extraction of cosmological pa- tQ make meaningful de-projections of gas temperature and den- 

rameters from large survey samples (X-ray, optical, or in the ra- shy profiles using SZE imaging data from multi . pixel bolometer 

dio/mm wavebands) relies critically on our understanding of dif- arraySi in combination with x . ray data . The A PEX-SZ exper- 

ferent mass observables, which depends on the detailed physical iment (Dobbfj et al 2QQ ^ Halverson et al 20 09) employs one 

processes affecting constituent gas and galaxies. of the firfjt such powerful mu lti-pixel Transition-Edge Sensor 

(TES) bolometer cameras, and a joint analysis of the ICM prop- 
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erties using SZE and X-ray data has been presented by Nord et 
al. (2009, Hereafter NBP09) for the massive cluster Abell 2163. 

In this paper we use the de-projection method used in 
NBP09 on the prototypical relaxed cluster Abell 2204. Our non- 
parametric analysis does not rely on any prior physical models in 
the construction of temperature and density profiles (apart from 
the assumption of spherical symmetry), hence the results are not 
based on parametric model fits. We measure the ICM entropy 
profile, as well as demonstrate the decrease of the ICM temper- 
ature in the cluster outskirts, first time from an SZE imaging 
data and independently from the X-ray spectroscopy. The de- 
rived ICM and cluster properties are compared with available 
X-ray and lensing results to highlight the level of accuracy of 
this independent method. 

Joint SZE/X-ray de-projection analysis is expected to be- 
come a standard tool in the near future for understanding the 
ICM physical state, as large numbers of resolved SZE maps 
will be available from the new generation SZE experiments. Our 
analysis assumes the gas to be in thermal equilibrium to model 
its physical properties, but presence of multi-phase ICM due to 
gas clumping will drive the electron temperature lower than the 
ion temperature in the electron-ion plasma (Evrard et al. 1996, 
Nagai et al. 2000). Recent hydro-simulations by Rudd & Nagai 
(2009) have shown, with a limited sample of halo models, that 
this deviation is small (about 5%) near rioa for a relaxed clus- 
ter. Joint SZE/X-ray analysis using interferometric measurement 
of the SZE with OVRO/BIMA (Reese et al. 2002) has already 
shown that clumping effects are not large in the cluster interior 
(within rsoo). Jia et al. (2008) have demonstrated the effect of the 
gas clumping on SZE and X-ray derived gas temperatures, and 
also found that these two quantities are in very good agreement 
within rsoo for the massive relaxed cluster RXC J2228. 6+2036. 
But at large radii the gas should get clumpier, due to the onset 
of filamentary structures. One vital goal for sensitive imaging of 
the SZE signal using wide-field, multi-pixel bolometer cameras, 
and its combination with the X-ray and weak-lensing measure- 
ments, will be to provide an ultimate tool for measuring the gas 
clumping and thermodynamic state near the cluster virial radius, 
to give a dynamic view on the growth of clusters through accre- 
tion. 

1.1. Previous SZE/X-ray joint modeling 

Due to the unavailability of resolved SZE images most of pre- 
vious SZE/X-ray joint analysis studies have been limited to an- 
alytical or numerically simulated cluster models with idealized 
noise properties. Zaroubi et al. (2001) considered a method for 
reconstructing the triaxial structure of clusters based on Fourier 
slice theorem and applied it to a set of cluster simulations. Lee 
& Suto (2004) also considered de-projection method combining 
SZE and X-ray data and applied to analytical cluster models. 
Puchwein & Bartelman (2006) have employed the Richardson- 
Lucy de-projection technique to reconstruct the ICM and probe 
the dynamical state of clusters from simulations, and Ameglio et 
al. (2007) used a joint SZE/X-ray likelihood function maximiza- 
tion using a Monte Carlo Markov Chain (MCMC) for a similar 
objective. 

Modeling ICM properties from real SZE observations has 
been limited mainly to isothermal /^-models (Cavaliere & Fusco- 
Femiano 1978). Holzapfel et al. (1997), Hughes & Birkinshaw 
(1998) used isothermal models to constrain the Hubble pa- 
rameter from observations of the clusters Abell 2163 and 
CL 0016+16, respectively, and later Reese et al. (2002) ex- 
tended this analysis to a sample of 18 clusters detected by 



OVRO/BIMA. De Filippis et al. (2005) used published SZE 
decrement values and X-ray imaging data to constrain the tri- 
axial structure of clusters using isothermal jfi-models. Zhang & 
Wu (2000) similarly used the /3-model to combine SZE and X- 
ray data to derive central gas temperature in clusters. A more 
detailed parametric modeling has been done by Mahdavi et al. 
(2007) for the cluster Abell 478, using simultaneous fits to the 
X-ray, lensing and SZE data assuming parametric models for 
dark matter, gas and stellar mass distribution, and hydrostatic 
equilibrium. 

Yoshikawa & Suto (1999) first used Abel's integral inver- 
sion technique, originally proposed by Silk & White (1978), 
for a non-parametric reconstruction of radial density and tem- 
perature profiles using analytical and simulated cluster models. 
More recently Yuan et al. (2008) has extended this method for 
the most X-ray luminous cluster RXC J 1347.5-1 145 using pub- 
lished /?-model fit values from SZE and X-ray measurements. 
Extrapolation of the density and temperature profiles to the clus- 
ter outskirts based on such parametric modeling can be problem- 
atic, in particular for clusters with a very peaked central emission 
such as RXC J1347.5-1 145. Additionally, this cluster is consid- 
ered to be a merging system (Cohen & Kneib 2002) where the 
assumptions of spherical symmetry and hydrostatic equilibrium 
may not be valid. The nearest approach to non-parametric mod- 
eling was made by Kitayama et al. (2004) for the same cluster, 
RXC J 1347.5-1 145, using a beta-model density profile to fit the 
X-ray surface brightness and obtaining fitted temperature values 
separately in each radial bin from their SZE imaging data. The 
small extent of their SZE map (less than 2 arcmin) limited the 
temperature modeling again to the cluster core region. 

1.2. Scope of the present work 

In this paper we apply the non-parametric ICM modeling based 
on Abel's integral inversion technique, as presented in NBP09, 
to the well studied and dynamically relaxed galaxy cluster Abell 
2204 (z = 0.1523, L x = 26.9 x 10 44 h~l erg s" 1 in the 0.1 - 2.4 
keV band, T x = 7.21 ± 0.25 keV; Reiprich & Bohringer 2002). 
The only assumptions in this analysis are spherical symmetry for 
reconstructing temperature and density profiles, and hydrostatic 
equilibrium (HSE) for reconstructing the total mass profile. The 
primary aim is to confirm the validity of this method for mod- 
eling the ICM distribution and cluster mass - and compare the 
results with those obtained from deep X-ray spectroscopic and 
weak lensing data - in a cluster where the assumptions of spher- 
ical symmetry and HSE are generally accepted to be valid. 

We compute the Chandra spectral temperature profile with 
the latest calibration updates and compare it with the SZE- 
derived temperature profile. In contrast to the X-ray spectro- 
scopic measurements from Chandra, the SZE-derived ICM tem- 
perature measurements near the cluster virial radius are con- 
strained primarily by the statistical uncertainties in the SZE data. 
This fact demonstrates the potential for stacking the SZE signal 
of several relaxed clusters to put tighter constraints on the slope 
of the gas temperature profile in the cluster outskirts (Basu et 
al., in preparation). For a single cluster (Abell 2204), the depth 
in the APEX-SZ map allows us to model the temperature profile 
with meaningful errors up to ~ 80% of the cluster virial radius 
(which we take to be r2oo, the radius within which the mean total 
density is 200 times the critical density). 

From density and temperature profiles we derive other physi- 
cal properties like total gravitational mass, gas mass fraction and 
the gas entropy index. The total mass modeling provides a quan- 
titative comparison with the published X-ray and lensing results. 
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The modeling of the gas entropy profile from SZE/X-ray imag- 
ing data is a first, and we compare the central entropy values 
of two clusters with different morphologies, A2204 and A2163 
(APEX-SZ analysis of the latter was presented in NBP09). This 
comparison shows how the gas entropy in the cluster core de- 
rived from SZE/X-ray joint modeling can be used to infer the 
dynamical state of clusters without the need for X-ray spec- 
troscopy. A further comparison of the baryonic fraction of the 
ICM between A2204 and two other dynamically complex clus- 
ters detected by APEX-SZ (Bullet and A2163) illustrates a sta- 
tistically significant difference of / gas inside ^on- 
All the scientific results in this paper are computed from the 
radial profiles of two observables: the SZE temperature decre- 
ment at 150 GHz, and the X-ray surface brightness in the 0.7 - 2 
keV band of XMM-Newton. In §2 and §3, we describe the map 
making and radial profile extraction steps from the X-ray and 
SZE data, and discuss the different systematic uncertainties as- 
sociated with each profile. §4 describes Abel's integral inversion 
method and presents our primary results in the form of the ra- 
dial density and temperature profiles. In §5 we present the other 
derived quantities like gas entropy and the total cluster mass pro- 
files, and list the conclusions in §6. 

We use the currently favored ACDM cosmology with the 
following parameters: Q m = 0.27, Qa = 1 - 0. m = 0.73, and 
the Hubble constant H Q = 70 km s _1 Mpc" 1 . At redshift of 
Z = 0.1523, the angular diameter distance of Abell 2204 is 541.6 
Mpc. To put the radial profiles in perspective using the charac- 
teristic cluster radii, we adopted the maximum likelihood NFW 
fit parameters from Corless et al. (2009), M 2 oo = 7.1 x 1O 14 M 
and c = 4.5, which gives r2oo = 1-76 Mpc (11.2'), ^500 = 1-16 
Mpc (7.3') and r 25 oo = 0.51 Mpc (3.2'). 

2. Extraction of the X-ray surface brightness profile 

This section describes the basic data analysis steps for X-ray 
map making, and the method for extracting the radial profile. A 
brief description of the analysis method is provided below, refer 
to references for further details. We discuss the main source of 
the X-ray systematic error caused by particle background, that is 
incorporated in the analysis. 

2.1. XMM-Newton observation and data reduction 

A2204 was observed by the XMM-Newton EPIC camera with 
medium filter in the full frame mode (ID: 0112230301). After 
carrying out the screening procedure (e.g. Zhang et al. 2008, 
hereafter ZF08) to filter flares, we obtained 17.5 ks, 18.5 ks 
and 14.3 ks clean exposure for the MOS1, MOS2 and pn in- 
struments. For pn data, the fraction of the out-of-time (OOT) 
events caused by read-out time delay is 6.30%, and a simu- 
lated OOT event file is created to statistically correct for this. 
The SAS command "edetect_chain" was used to detect point- 
like sources, which were subtracted before further data reduc- 
tion. The vignetting correction to the effective area is accounted 
for by the weight column in the event lists. Geometric factors 
such as bad pixel and gap corrections are accounted for in the 
exposure maps. We choose the XMM-Newton blank sky accumu- 
lations in the Chandra Deep Field South (CDFS) as background. 
The background observations were processed in the same way 
as the cluster observations. The CDFS observations used the thin 
filter, while the A2204 observations used the medium filter. The 
background of the A2204 observations is thus different from the 
CDFS using the thin filter at energies below 0.7 keV. Therefore 
we performed all the analysis at energies above 0.7 keV, in which 




Fig.l. 10 x 10 arcmin XMM-Newton MOS1 image of A2204, 
flat-fielded and smoothed with a 12" wide gaussian kernel. The 
overlaid contours (white) are from the APEX-SZ measurement; 
contour steps are -2, -4, -6 and -8 cr, and the SZ image resolution 
is 1 arcmin. The X-ray image is made in the 0.7-2 keV band, 
binned in 4" pixels. The color bar is in a logarithmic scale of 
0.1- 400 counts per pixel. The black cross in the center denotes 
the flux weighted X-ray center used for the surface brightness 
profile extraction, and is within 4 arcsec of the SZE peak location 
obtained from spherical /J-model fit. 



the difference of the background is negligible. The image of 
A2204 is shown in FigQ] 

2.2. X-ray profile extraction 

The 0.7-2 keV band is used to derive the surface brightness 
profiles. This ensures an almost temperature-independent X-ray 
emission coefficient over the expected temperature range. The 
width of the radial bins is 2" . An azimuthally averaged surface 
brightness profile of the CDFS is derived in the same detector 
coordinates as for the target. The count rate ratios of the tar- 
get and CDFS in the 10-12 keV band and 12-14 keV band for 
MOS and pn, respectively, are used to scale the CDFS surface 
brightness. The residual background in each annulus of the sur- 
face brightness is the count rate in the 0.7-2 keV band of the 
area scaled residual spectrum obtained in the spectral analysis. 
Both the scaled CDFS surface brightness profile and the resid- 
ual background are subtracted from the target surface brightness 
profile. The background subtracted and vignetting corrected sur- 
face brightness profiles for three detectors are added into a single 
profile, and re-binned to reach a significance level of at least 3- 
cr in each annulus out to r < 9 arcmin. The particle-induced 
background varies by less than 10% comparing the background 
observations. Therefore the dispersion of the re-normalization 
of the background observations is typically 10%. We take into 
account a 10% uncertainty of the scaled CDFS background and 
residual background. The resulting profile is shown in Fig|2] and 
values are given in Table 1 . 

Note that raising the upper cut of the energy band does not 
provide dramatic improvement in the signal-to-noise ratio in the 
surface brightness profile in the cluster outskirts (> 3'), where 
the gas has a lower temperature and therefore does not contribute 
significant X-ray photons at high energies. 
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Fig. 2. X-ray surface brightness profile from the XMM-Newton 
map in black squares, converted to physical units. The profile 
has been convolved and re-binned to match the APEX-SZ res- 
olution. The errors include the systematic uncertainties due to 
background modeling and are incorporated in the de-projection 
analysis. Over-plotted in green are the surface brightness values 
prior to re-binning. 



The X-ray surface brightness profile obtained from the above 
procedure is convolved with a one-arcmin gaussian kernel to 
bring its resolution to the same level as for the APEX-SZ raw 
image. This smoothing raises the S/N ratio, particularly in the 
cluster outskirts. Additional re-binning is performed to conform 
the X-ray profile with the SZE data, since the latter is averaged 
in wide annular bins in the cluster outskirts to keep the statis- 
tical uncertainties under control. The widths of the radial bins 
are given in Table 1, where the central bin corresponds to the 
central 1 arcmin circle in each map. The resulting X-ray surface 
brightness profile after re-binning is shown in Figf2] (squares). 
The errors indicate 1 <x uncertainty from the added poisson er- 
rors of the target and CDF surface brightness profiles, plus sys- 
tematic uncertainties due to the variation of the particle-induced 
background. Beyond 6 arcmin radius the systematic uncertainty 
starts to dominate, but we still have a S/N of > 3 within 9 arcmin 
radius. As we will discuss in the next section, our results are 
currently dominated by the statistical uncertainties in the SZE 
imaging at this radius, therefore we are not limited by X-ray 
systematics in the modeling of one single cluster. 

2.3. X-ray imaging vs. spectral spatial ranges 

Here we briefly highlight the advantage of the SZE/X-ray de- 
projection method to obtain the ICM temperature out to large 
radii, in comparison with the X-ray spectral measurements (in 
particular from XMM-Newton and Chandra). To constrain the 
gas temperature to an uncertainty smaller than 10% from the 
X-ray spectra for such a hot cluster, one needs typically S/N 
> 150 in the 0.7-7.8 keV energy band (e.g. Zhang et al. 2009). 
Therefore such temperature measurements are typically limited 
to the central regions of clusters (up to a radius between r25oo 
and rsoo, less than half of the cluster virial radius, see ZF08). 
Recently data from the Suzaku satellite have been used to mea- 
sure gas temperature beyond rgoo for a few clusters (Fujita et al. 
2008, Reiprich et al. 2009, George et al. 2009). However, these 
observations are expensive and limited to a few nearby (z < 0.2) 
clusters only. SZE/X-ray joint modeling can overcome this is- 
sue, by using X-ray surface brightness to provide primarily a 
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Fig. 3. 10 x 10 arcmin APEX-SZ map of Abell 2204, overlaid 
with XMM-Newton X-ray contours in steps of 3cr in logarithmic 
scale and smoothed to the APEX-SZ resolution. The SZE map 
has been deconvolved to the beam scale to reconstruct the full 
signal (see text). The FWHM of the APEX-SZ beam is shown in 
the lower left. 



constraint on the gas density, and then obtaining the temperature 
from SZE data. This easily allows for measuring the gas tem- 
perature at the outer radii where the X-ray S/N is low, e.g. 3-5. 
Thus ICM modeling up to the cluster virial radius can be done, if 
the systematic uncertainties in both the X-ray and SZE imaging 
are controlled, and the SZE statistical uncertainties are brought 
down. 



3. Extraction of the SZE temperature decrement 
profile 

This section describes the basic reduction and map making steps 
for the APEX-SZ data. The analysis is very similar to that of 
NBP09, which should be consulted for further details. Here we 
emphasize the construction of a set of SZE temperature decre- 
ment profiles, all consistent with our APEX-SZ measurement, 
that we use to estimate the uncertainties in the de-projection 
analysis. A similar approach was also used in NBP09, but the 
details of SZE profile construction and de-projection procedure 
were not discussed. 

3.1. APEX-SZ observation and map making 

Abell 2204 was observed with the APEX-SZ camera in May 
2008 and April 2009, with roughly 80% of the observing time 
spent in the 2008 run. The usable data on the target amounts to 
approximately 10 hours, divided between scans of 20 minutes 
duration. The primary calibration source at 150 GHz was Mars, 
and secondary calibrators were Neptune and RCW38. Details of 
the observing technique and data calibration for APEX-SZ are 
given in Halverson et al. (2009, hereafter HL09). We also re- 
fer to the Fig.l of that paper for an illustration of the scanning 
pattern. 

The reduction process is optimized for the circular drift scans 
employed for the observation. After eliminating detectors with 
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low optical response, correlated atmospheric noise is removed 
by subtracting the median signal across the good channels of 
the array at each time step after a temporary normalization step 
has been carried out. Additional reduction steps such as despik- 
ing and de-glitching are used, but affect only a small amount of 
data. Circular subscans are baselined, by defining the subscans 
consisting of 3 full circles and then applying a fifth order poly- 
nomial, which corresponds to a low-pass spatial filter (affecting 
spatial scales only marginally larger than those filtered out by the 
circular scan pattern itself). For each scan, a map with 15" x 15" 
pixels is constructed, weighting the data by the inverse rms at 
the position of each pixel in each scan. The result we refer to 
as the "raw map", and the radial profile made from this map is 
shown in Fig]4] In parallel, a bright point source convolved with 
the instrument beam (obtained from fitting the Mars scans out 
to a 4.5' radius) is processed by an identical pipeline to obtain 
the transfer function (see HL09), which is used to perform the 
deconvolution. 

The deconvolution of the map is performed iteratively in map 
space as described by NBP09 and discussed in more detail by 
Nord (2009, PhD thesis). The process essentially reconstructs 
the cluster signal as the sum of many point sources as seen by the 
instrument beam. The final deconvolved map is shown in Figj3] 
overlaid with the X-ray surface brightness contours. The noise 
on scales equal to the APEX-SZ beam is 44 /uKcmb in the central 
region of the map, corresponding to a peak signal-to-noise ratio 
of 8.5. 

The outer contours of the APEX-SZ map with low signal- 
to-noise ratio shows an elliptical shape. This is most likely the 
result of unfiltered noise on scales of several arcminutes (but see 
Corless et al. 2009 for a discussion on the triaxial dark matter 
halo in this cluster). We perform spherical and elliptical isother- 
mal j6-model fits to the SZE map, which yield identical val- 
ues for the SZE emission center; (RA,Dec) = (248.196,5.577). 
These coordinates are within 4 arcsec of the flux-weighted X- 
ray emission center, which is defined iteratively though a series 
of concentric circles in the X-ray map (see Zhang et al. 2010, 
§2.3). This 4 arcsec offset is comparable to the pointing accu- 
racy of the APEX telescope at 150 GHz. This provides additional 
confirmation for the relaxed morphology of this cluster, to apply 
spherical de-projection using a common SZE/X-ray center. 

3.2. SZE profile extraction and noise properties 

To estimate how uncertainties in the SZE map are propagated 
through our analysis, we compute a set of (typically 100) decon- 
volved SZE profiles by applying the above map making process 
on different noise realizations. The resulting profiles are used to 
compute all the relevant cluster properties (i.e. profiles of gas 
density and temperature, and thereafter mass and entropy pro- 
files). 

We first obtain a set of "jack-knived" noise maps, by ran- 
domly selecting half of the scan maps and inverting their signs, 
and then co-adding them with the remaining scans (see HL09 
for details of this technique applied to APEX-SZ data). This re- 
moves all astrophysical signals but retains the noise structures 
unfiltered by the pipeline. A random realization of the Cosmic 
Microwave Background (CMB) sky is added to these noise maps 
to estimate the contamination of the SZE signal by the CMB 
(§3.3). An azimuthally symmetric cluster map is made from the 
radial profile of the raw map, and added to the jack-knived noise 
maps. The thus simulated raw maps are deconvolved using the 
transfer function, and a set of radial SZE profiles are obtained 
from the final maps. The scatter in this set of profiles constitutes 



0.1 F 




-0.4 r / i 

-0.5 F', , , i , , , i , , , i , , , i , , , i , , , i , 
2 4 6 8 10 12 

Angular distance (arcmin) 

Fig. 4. Radial profile of the SZE temperature decrement at 150 
GHz. The profile obtained from the raw reduction is shown by 
the blue diamonds, and the mean deconvolved profile by the 
red squares. The errors on the deconvolved profile represent the 
l<r scatter in the set of profiles used in the de -projection anal- 
ysis (see text). These errors are correlated especially within the 
narrow central bins. This reduces the significance of the peak 
decrement as seen directly from the deconvolved map (Pigf3j- 
We show a random set of 5 deconvolved profiles (truncated at 8 
arcmin) to illustrate the noise correlation. 



the total statistical uncertainty in the SZE measurement, shown 
in Figg] 

The primary motivation for constructing a set of SZE radial 
profiles from jack-knived noise maps is to incorporate the ef- 
fect of noise correlation occurring due to the presence of unfil- 
tered noise structures in the map, which typically have scales 
much larger than the APEX-SZ beam. This reduces the signifi- 
cance of the detection of the SZE signal, in particular if narrow 
binning is used. This is clearly seen from the errors on the fi- 
nal deconvolved radial profile in FiglU which are about 50% 
larger than the errors computed from the variance in each an- 
nular bin in the deconvolved cluster image (Figf3j. The total 
statistical uncertainties are then easily propagated through the 
de-projection analysis by computing the relevant physical quan- 
tities for each profile and measuring their scatter in each radial 
bin. This method also makes sure that the numerical errors com- 
ing from the de-projection method are not artificially enhanced 
(see details in §4.3). 

Additionally, the choice of 1 arcmin binning in the central 
region of the SZE map leads to a correlation between the adja- 
cent bins due to PSF smearing. The measured solid angle of the 
APEX-SZ beam is 1.5 arcmin 2 , and 30% of the beam power is 
in the near sidelobes outside the best fit gaussian beam with full- 
width at half maximum (FWHM) of 58" (Reichardt et al. 2009). 
Density and temperature de-projection based on a proper PSF 
deconvolution is not attempted in this paper; we simply note that 
the errors in the narrow central bins are possibly under-estimated 
by a small amount, with an overall downward bias in the mea- 
sured gas temperature in these bins. The current choice of nar- 
row binning inside r2soo of the cluster is motivated by the aim 
of demonstrating the compatibility of our temperature and mass 
profiles with published X-ray results. It is also not desirable to 
smooth out the effect of the central cool core of Abell 2204. 

The amount of correlation present between different radial 
bins is easily computed by means of the correlation matrix. We 
compute the correlation matrix for the set of 100 deconvolved 
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SZE profiles and find that the 4 central bins are almost fully 
correlated (p > 0.8), while the outermost bins have little cor- 
relation. For uniform binning of the SZE profile (12 bins of 1 ar- 
cmin each) there is correlation between adjacent bins all across 
the profile due to large scale noise structures, and in addition the 
bins near the center are correlated more strongly than the oth- 
ers due to PSF smearing. However, uniform binning is not used 
at large radii for extracting information out to the very low S/N 
regions of the SZE map, in a cluster whose detection signifi- 
cance is lower than those presented previously from APEX-SZ 
(HL09, NBP09). When fewer broad bins are used, the correla- 
tion becomes negligible as can be expected (p < 0.1 for 4 equal 
bins), but this is not used either as we are interested in the clus- 
ter cool core. This correlation pattern for any radial binning will 
propagate through all the other derived cluster quantities (like 
temperature, total mass and entropy bin values). Additionally, 
the noise will not go down as expected when averaging several 
bins due to correlated errors, therefore we re-bin the original set 
of deconvolved maps for computing errors on averaged values. 

3.3. Sources of systematic errors 

The deconvolution method used in making the final SZE map 
can introduce systematic bias in the final profiles. As in NBP09, 
a series of simulations is performed by adding artificial cluster 
models (^-profiles) to jack-knived noise maps and passing them 
through the reduction pipeline. These are then deconvolved us- 
ing the transfer function, and the resulting profiles are compared 
with the input y6-models. The effect is a systematic lowering of 
the cluster signal at large radii due to flux loss, by as much as 
40% at f2oo (this number is true only if the real cluster profile 
follows an isothermal /3-model). This error is considerably lower 
than the intrinsic statistical uncertainties on the profile, which in 
case of A2204 is almost 100% at r200- Resulting systematic un- 
certainties on the ICM temperature profile are discussed in §4.3, 
after describing the de-projection method. 

For a large cluster like A2204 (virial radius ~ 12 arcmin), 
the temperature anisotropies in the CMB are a major source of 
confusion. Following the same prescription as in NBP09, we 
attempt to quantify this by making multiple realizations of the 
CMB sky using the HEALpix software (Gorski et al. 2005), and 
adding these to the jack-knived noise maps before performing 
noise simulations. The additional scatter in the resulting radial 
profiles is 14% at r 5 oo and roughly twice as large at r2oo, again 
less than the statistical errors in the APEX-SZ measurements. 
The systematic uncertainty on the APEX-SZ measurement aris- 
ing from calibration errors is of the order of 5% (HL09). 

Other systematic errors in the SZ map can arise from un- 
resolved point sources (radio or sub-millimeter galaxies) and 
galactic dust emission, which we have ignored. There is no indi- 
cation of any point like sources in the 150 GHz SZE map. The 
NVSS radio catalog (Condon et al. 1998) lists a 70 mJy radio 
source at 1 .4 GHz approximately 9" from the cluster X-ray cen- 
ter. After subtracting the best fit /3-model from the raw map the 
rms noise at the map center is 2.2 mJy/beam, and no indication 
of a point source is seen in the residual raw map. The IR lu- 
minosity of the central brightest cluster galaxy (BCG) in Abell 
2204 is reported by Quillen et al. (2008). The corresponding dust 
thermal emission at 150 GHz (< 0.1 mJy/beam) is much be- 
low the noise level at the map center, and the downward bias 
in the measured gas temperature at the cluster center can be 
ignored. Knudsen et al. (2008) found a bright sub-millimeter 
galaxy, SMM J163244.7+053452, in the field of A2204 at a 
distance of 39 arcsec N-W from the cluster center. Its 850 //m 



flux density is 22.2 ± 4.9 mJy, with estimated magnification of 
li = 3.4. Assuming a spectral index a = 3, where S v oc v a , 
this source will produce a flux density of roughly 1.7 mJy at 150 
GHz, corresponding to a temperature increment of 34 fiK for the 
APEX-SZ beam. This is lower than the noise rms at the position 
of this galaxy in the map. 

4. De-projection of radial density and temperatures 

The three-dimensional (de-projected) density and temperature 
profiles are obtained directly using Abel's integral inversion 
method (as in NBP09), with the assumption of spherical sym- 
metry. Although proposed nearly three decades ago for joint 
SZE/X-ray analysis (Silk & White 1978), this method has re- 
mained largely unused. One possible reason for this limited 
application might be due to its numerical instability, as it in- 
volves computing derivatives at each point on the observed pro- 
files. We have utilized the noise correlation in the real SZE data 
to partially overcome this problem, which makes Abel's inver- 
sion technique a particularly simple and intuitive method for 
de-projection. Unlike the standard "onion-skin" method of de- 
projection used in X-ray spectral analyses (Kriss et al. 1983), 
Abel's inversion is not dependent on the choice of the outer- 
most bin. The strong anti-correlation in the de-projected temper- 
ature values between adjacent bins, a numerical artifact found in 
several geometrical de-projection techniques (see Ameglio et al. 
2007), is also not significant. 

4. 1 . Method for de-projection 

For the de-projection analysis, the SZE temperature decrement 
can be written as the integral of the electron pressure along the 
line of sight as 

AT(R) = 2A SZE f g{x,T e )n e (r)T e {r) ^ (1) 
Jr Vr 2 - R 2 

where Asze = Tqmb (fco" T /m e c 2 ), r is the physical radius from 
the cluster center, R = D^O where 8 is the projected angular dis- 
tance on the sky and D A is the angular diameter distance. T e (r) 
and n e (r) are the electron gas temperature and density radial 
profiles. g(x, T e ) is the frequency dependence of the SZ signal, 
in which the gas temperature dependent relativistic correction 
terms have a small contribution (~ 5% at 150 GHz for a 10 keV 
cluster). Therefore, we can neglect the radial temperature de- 
pendence in g(x, T e ) and incorporate a fixed temperature value 
g(x, r e =8.26 keV) into the Asze factor, following the mean X- 
ray temperature from Arnaud et al. (2005). Note that we are ig- 
noring any contribution from the kinematic Sunyaev-Zel'dovich 
effect (kSZE), as its contribution is likely to be much less than 
our imaging uncertainty. 

In a similar way, the X-ray surface brightness profile can be 
written as 

S * ( V = TTi * n 2 e (r) A H (T e (r)) = ■ (2) 

4tt(1 +z) 4 J r Vr 2 - R 2 

We compute the value of the X-ray emissivity function A H (T e ) 
in each radial bin with the MEKAL code in XSPEC (Mewe et al. 
1982, Kaastra 1992), using models for metallicity and tempera- 
ture radial profiles obtained from the spectral measurements of 
ZF08. The actual measured metallicity values within 0-3 arcmin 
radius changes from 0.5 Z Q to 0.3 Z G , corresponding change in 
A#(r) is 12%. The weak temperature dependence of the soft X- 
ray emission in the 0.7 - 2 keV energy band makes our results 
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practically insensitive to any temperature model used in calcu- 
lating A H (T e ), as noted in NBP09. 

Using Abel's integral equation, equations (Q]) and (O can be 
inverted to obtain joint radial density and temperature profiles 
(Yoshikawa & Suto 1999) 
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Equations (0) and d?} are integrated numerically by summing 
in radial bins from i m { n to ; max , where / max is the index for the 
outermost bin, and / m ; n corresponds to r/D^. 

To show that our analysis results do not depend on any a pri- 
ori knowledge of the radial temperature profile, we tried two al- 
ternative approaches for the computation of the emissivity func- 
tion in addition to the MEKAL model. We used a mean value of 
the X-ray temperature in all bins to compute A//(r c (r) = Tx), 
where Tx = 8.26 keV (Arnaud et al. 2005). Alternatively, we 
used a weak power-law dependence of the emissivity function 
on the gas temperature, as Afj(T e (r)) oc T e (ry 1 ^ 6 . The second 
assumption gives excellent approximation to the X-ray emissiv- 
ity values from the MEKAL code if we assume the bulk of the 
cluster gas has a temperature in the range 2-14 keV. The de- 
projected ICM density and temperature profiles from all three 
methods, after combining with the SZE radial profiles, are es- 
sentially identical given the statistical error in each radial bin. 
This confirms the fact that the use of the soft-band X-ray data in 
our analysis is primarily providing the constraints on gas density 
profile, whereas the gas temperature constraints come from the 
SZE measurement. 



4.2. Radial profiles for gas density and temperature 

The results for de-projection of density and temperature profiles 
for A2204 are shown in Fig|5] and the corresponding values with 
their errors are given in Table 1 . Note that the uncertainties on 
the X-ray surface brightness profile due to the background mod- 
eling are included in the results, but the additional (small) sys- 
tematic uncertainties from X-ray flux calibration are neglected, 
which likely produces an under-estimation of the errors on the 
density values in the inner bins. In the upper panel of Fig [5] we 
overplot the density profile obtained by ZF08 by fitting a double 
/?-model to the X-ray surface brightness. This density profile is 
XMM-Newton PSF corrected, and the common B slope parame- 
ter is obtained by fitting the outer component (see A. 3 in ZF08). 
The rms fractional errors shown below Fig [3] are computed as 
X = (bin density - model density) /bin error. Except for the inner 
arcminute where the X-ray brightness profile is extremely cuspy, 
the double beta model provides a good fit to our de-projected 
bin densities. This follows from the fact that in the 0.7-2 keV 
energy band the X-ray surface brightness is practically indepen- 
dent of the gas temperature. A similar argument had been used 
by Kitayama et al. (2004) while modeling the gas density profile 
with the X-ray derived y6-model to obtain best fit radial tempera- 
ture values in RXC J 1347-1 145. 

In the lower panel of Fig|5]the radial temperature profile is 
shown. There is a clear indication of the cluster cool core from 
APEX-SZ data; the temperature drops almost by a factor 3 from 
500 kpc to 100 kpc radius. This is in contrast with the temper- 
ature profile for A2163 in NBP09, which could be fitted with a 
single isothermal profile at all radii within the lcr uncertainties 
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Fig. 5. Top panel: The de-projected electron density with cor- 
responding errors. The dashed line is the density profile from 
ZF08, obtained from fitting a double y6-model to the X-ray sur- 
face brightness. The rms fractional differences from the yS-model 
fit are shown below. Bottom panel: The de-projected temperature 
values and their lcr statistical errors from the SZE measurement. 
Over-plotted data points in red diamonds are from XMM-Newton 
analysis by ZF08 (see §4.4). The arrows in both plots mark the 
characteristic cluster radii, and the horizontal "error bars" the 
width of the bins. 



of the bin values. There is also a strong indication of a decreas- 
ing ICM temperature beyond its peak at ~ 500 kpc. However, 
the temperature value at the last radial bin at r2oo is essentially 
an upper limit, there is no statistically significant SZE signal at 
this radius given the current noise level in the map. 

The decreasing temperature profile in the cluster outskirts 
can be made clearer by re-binning the APEX-SZ data. As noted 
in §3.2, the errors in the bin values are correlated and do not 
average down as expected in random gaussian noise. Therefore, 
we re-bin the original set of deconvolved maps to compute the 
bin errors. The result is shown in Figj6] where we have divided 
the data in only two bins, excluding the central 3' of the map. 
A decrease in gas temperature from its peak value is supported 
at 98% confidence level (2.3<x). To put this temperature slope in 
perspective, we overplot in Fig|6]the results from recent X-ray 
observations and numerical simulations of clusters, scaled to the 
values for A2204. The solid line is the mean spectroscopic tem- 
perature profile in cooling core clusters (Vikhlinin et al. 2005), 
and the gray shaded region is the average profile of the cooling 
core clusters from ASCA with their lcr dispersion (Markevitch 
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Fig. 6. The APEX-SZ measurement of the outer temperature 
profile in the cluster A2204, as compared to simulations and X- 
ray spectral measurements. The two data points represent the 
SZE-derived temperature values with lcr errors; a decrease in 
gas temperature is supported at 98% confidence level. The red 
dashed region represents a fit using UTP profiles from cluster 
simulations by Hallman et al. (2007), with lcr errors. The gray 
shaded area shows the average profile of cooling core clusters 
from ASCA (Markevitch et al. 1998), and the blue solid line is 
the Chandra observation of cooling core systems (Vikhlinin et 
al. 2005). 



et al. 1998). The SZE radial temperature is statistically consis- 
tent with both these measurements, although it appears to in- 
dicate a steeper slope. The SZE-derived temperature slope also 
appears steeper than the Universal Temperature Profile (UTP) 
fit from numerical simulations of relaxed cluster (Hallman et al. 
2007), shown in the dashed line with the hatched region for the 
lcr uncertainties in the UTP fit values. Again, SZE measurement 
from one cluster is not yet adequate to provide a quantitative 
comparison with the numerical simulations, but a stacking anal- 
ysis of several relaxed clusters can be expected to yield a mean- 
ingful comparison by lowering the statistical noise. 

As mentioned in the introduction, the two quantities T e and 
r g as are used synonymously in this paper, where the latter is 
defined as r gas = (n e T e + njTi)/(n e + «,). They can differ if 
the post-shock equilibrium timescale between the electrons and 
ions is large, especially in the low density environment near 
rioo (where n e ~ 10~ 4 - 10~ 5 cm -3 ). Recently, Rudd & Nagai 
(2009) have provided quantitative estimate for this bias from 
cosmological hydrodynamic simulations of galaxy clusters, and 
found that for relaxed clusters (CL 104 in their simulations, with 
Tgas = 5.4 keV) T e can under-estimate the r gas by about 5% at 
r2oo- Therefore this effect can be ignored for our current analy- 
sis of a single cluster. Their results most likely present the upper 
limit of this bias, since non-adiabatic heating due to shocks in 
the cluster outskirts is neglected in their models. 



4.3. Errors on the de-projected profiles 

As seen from Eqns.(f3]l and the de-projection of density and 
temperature requires taking derivatives at each radial bin, which 
are the major source of introducing modeling errors onto the de- 
projected profiles. This fact may possibly have limited the appli- 
cation of Abel's inversion in the SZE simulations, using realis- 
tic mock observations with white noise. Although the high S/N 
imaging data from APEX-SZ with 1 arcmin resolution makes the 



application of Abel's inversion method feasible for the first time, 
propagating errors through a Monte-Carlo method will lead to 
a high and unphysical error level if the noise correlations be- 
tween radial bins are ignored. As noted by Yoshikawa & Suto 
(1999), pre-smoothing the data will reduce this error, but due to 
its model dependent nature we refrain from smoothing. It is also 
difficult to determine the degree of smoothing to be applied: a 
small smoothing kernel is insufficient to offset the numerical er- 
ror (in particular for the narrow central bins), whereas smoothing 
over several bins will make their errors artificially low and intro- 
duce modeling bias. 

The construction of a set of radial SZE profiles from jack- 
knived noise maps, described in §3.2, is used to overcome this 
problem. The X-ray systematic error due to uncertainties in the 
background modeling is treated as an amplified white noise, and 
random realizations of X-ray brightness profiles are obtained. 
Each of these X-ray profiles are then combined with one de- 
convolved SZE profile, and the de-projected density and tem- 
perature profiles are obtained. The scatter in each SZE profile is 
reduced by noise correlation, which keeps the numerical errors 
coming from Abel's inversion method at a minimum. Apart from 
density and temperature, profiles for all other cluster properties 
(like total mass, entropy) are obtained similarly: the scatter of 
the profiles measures the statistical uncertainties in each bin. The 
treatment of X-ray systematics as random noise is justified as the 
uncertainties in the current de-projected temperature values orig- 
inate almost entirely from the SZE measurement. For compari- 
son, estimating errors from a "blind" Monte-Carlo method treat- 
ing the SZE decrement value in each bin as independent gives 
temperature profile errors that are on average 2-4 times higher, 
thus making a demonstration of the decreasing gas temperature 
in the cluster outskirts impossible. 

The effect of SZE systematic errors on the gas temperature 
measurements are computed by methods described in §3.3. The 
relative amplitude with respect to statistical uncertainties and 
the radial dependence of the SZE systematic errors are simi- 
lar to those found for Abell 2163 in NBP09. That work pre- 
sented tabulated uncertainty values on both T e and n e . We ignore 
systematic uncertainties on gas density as it is much more ro- 
bustly constrained than the gas temperature. The systematic un- 
certainties on T e at r5oo due to confusion with the primary CMB 
anisotropies is +13%, and at r2oo it increases to nearly twice that 
amount. Irrecoverable loss of the SZE signal occurs for scales 
larger than the path of a single bolometer, and hence can not 
be recovered by the transfer function (see §5.2 in NBP09). This 
bias amounts to a systematic uncertainty of +15% at rsnn, and 
roughly +40% near r2oo (these numbers are accurate only in the 
context of an isothermal /3-model). For comparison, the statisti- 
cal uncertainty in the SZE map at r^a is about 45%, increasing 
to ~ 100% at r2oo- In Fig(7]we have shown the effect of adding 
systematic uncertainties (by adding in quadrature with statistical 
errors) on the measured gas temperature values. 

4.4. Comparison with X-ray spectral analysis 

In Fig|5]we also showed the de-projected X-ray spectroscopic 
temperature measurements for A2204 (ZF08), to provide a di- 
rect comparison between our results and those derived from X- 
ray analysis. There is partial overlap between these two pro- 
files within their lcr uncertainties, however, near r25oo the SZE- 
derived temperature is systematically higher. It is beyond the 
scope of this paper to present a detailed discussion on X-ray 
spectral analysis and its biases, but we mention the fact that 
for multi-temperature ICM in hot clusters the spectral analysis 



Basu et al.: Non-parametric ICM modeling from APEX-SZ imaging data 



9 



Table 1. The input X-ray surface brightness and SZE temperature decrement values at each radial bin, and the different de-projected 
quantities derived from these two observables. The errors shown are the l<x statistical uncertainties, including the systematic errors 
from the X-ray background modeling (the latter is included in the X-ray surface brightness errors). The first bin represents the 
central arcmin of the maps, corresponding to the APEX-SZ beam FWHM. 



Bin Radius 


Sx 


Ar sz 


n e {r) 


T e {r) 


Entropy^ 
(keV cm") 


M tot (< r) 


(arcmin) 


(t) 


(mK) 


(1(T 3 cm" 3 ) 


(keV) 


(10 14 A/ Q ) 



1 


0- 


-0.5 


10.8 ± 0.2 


-0.36 ± 


0.06 


24.4 i 


= 0.2 


2.7 ± 


0.5 


32 ±6 


0.03 ± 0.01 


2 


0.5 


- 1.5 


1.03 ±0.02 


-0.28 ± 


0.05 


5.30 ± 


0.04 


7.2 ± 


0.9 


240 ± 30 


0.54 ± 0.04 


3 


1.5 


-2.5 


0.18 ±0.007 


-0.18 ± 


0.04 


1.79 ± 


0.03 


10.4=1 


= 1.7 


710 ± 120 


1.72 ±0.12 


4 


2.5 


- 3.6 


0.06 ± 0.003 


-0.13 ± 


0.03 


0.86 ± 


0.02 


12.9=1 


= 2.2 


1430 ± 270 


2.96 ± 0.41 


5 


3.6 


-4.9 


0.02 ± 0.002 


-0.06 ± 


0.02 


0.47 ± 


0.02 


9.0 ± 


2.1 


1490 ± 370 


3.82 ±0.83 


6 


4.9 


-6.7 


0.01 ± 0.001 


-0.03 ± 


0.01 


0.27 ± 


0.02 


5.8 ± 


1.9 


1390 ±510 


3.64 ±0.61 


7 


6.7 


-9.2 


(2.8 ± 1.0) x 10" 3 


-0.014 ± 


0.007 


0.11 ± 


0.02 


4.0 ± 


2.5 


1700 ± 1200 


3.9 ± 1.3 


8 


9.2- 


- 12.8 


(1.0 ± 0.8) x 10" 3 


-0.010 ± 


0.006 


0.06 ± 


0.02 


5.7 ± 


5.3 


< 7500 c 


< 13.1^ 



" entropy index, defined as K = k^T e n\ 
b in units of 10~ 12 erg s _I cm -2 arcmin 
c upper limit at 68% confidence level 
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method can significantly underestimate (by up to 40%) the mass- 
weighted gas temperature (Mazzotta et al. 2004), and this effect 
is expected to be stronger near cluster cool cores where the line 
of sight crosses many temperature components. The low temper- 
ature value in the innermost bin from our measurement may be 
partially caused by the numerical uncertainty of taking deriva- 
tives at the inner edge of the profile, or APEX-SZ beam smear- 
ing. Snowden et al. (2008) considered the effect of XMM-Newton 
PSF smearing in analysis of this cluster, and gave a higher value 
of X-ray spectroscopic temperature near 1' radius. 

In order to avoid added complexities from the X-ray spec- 
tral de-projection, a simpler way is to make a projected (i.e. two 
dimensional) temperature profile from our measurements using 
an appropriate weighting scheme. The mean weighted value of 
the gas temperature along the line of sight can be computed 
as r pl qj = J WTdV I J WdV, where T is the de-projected gas 
temperature and W is the weight function. We use two dif- 
ferent weighting schemes: the standard emission weight with 
W = n 2 A(T) (using A(r) oc r" 1/6 as discussed earlier), and the 
weighting for a "spectroscopic-like" temperature as discussed by 
Mazzotta et al. (2004), using W = n 2 T ~ 3/4 . As seen in NBP09, 
the results are almost identical for these two methods, and pro- 
jection results only for the Mazzotta model are used for com- 
parison with the X-ray data. As can be expected, the effect of 
projection on the radial temperature profile is small when com- 
pared to the current statistical errors. 

For an accurate measurement of the X-ray spectroscopic 
temperatures, we have re-analyzed two Chandra observations of 
A2204 (Obs IDs 6104 and 7940), resulting in a total exposure of 
88 ks. While a temperature profile using these data has already 
been published (Sanders et al. 2009), the calibration update re- 
leased recently (CALDB 4) was expected to have a significant 
effect for this hot cluster. Therefore, the Chandra data was re- 
analyzed in the same way as described in Hudson et al. (2010); 
in addition, a correction for a possible difference in the cosmic 
X-ray background between source and blank sky observations 
was applied since we are also interested in low surface bright- 
ness cluster outskirts. The analysis with the new calibration re- 
sults in approximately 15-20% lower temperatures in the hot 
cluster regions (at < 7"25oo) as compared to Sanders et al. (2009). 

The results from the updated Chandra spectral analysis and 
the projection of the SZE/X-ray 3D temperature profile are 
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Fig. 7. Comparison of the projected gas temperature profile in 
A2204 deduced from APEX-SZ data (blue, solid boundaries) 
with X-ray spectral measurements from Chandra data (red, 
dashed boundaries). The Chandra spectral analysis has been 
re-performed taking the latest calibration update into account. 
The hatched regions show the Icr statistical uncertainties in each 
measurement. On top of that we overplot the total uncertainties 
in each method combining statistical and systematic errors in 
quadrature. 



shown in Figj7] The blue hatched region marks the ±lcr sta- 
tistical uncertainties around the mean SZE-derived temperature, 
and similarly the red hatched region shows the statistical uncer- 
tainties in the X-ray spectral analysis. Both results are in excel- 
lent agreement within their mutual uncertainties. But the point to 
note in Fig|7]is the relative increase in the statistical and system- 
atic errors in the Chandra and APEX-SZ measurements of the 
gas temperature. The statistical and systematic errors are added 
in quadrature and the total uncertainties are shown on top of 
the statistical uncertainties (white bordered regions). At r^m the 
Chandra measurement is already dominated by systematic un- 
certainties due to the background modeling; beyond that radius 
it is impossible to put meaningful constraints on the gas tem- 
perature using the current Chandra data. In contrast, the uncer- 
tainties on the SZE-derived temperatures are dominated by the 



10 



Basu et al.: Non-parametric ICM modeling from APEX-SZ imaging data 



statistical errors also at r2oo- At rsoo the ratio of statistical and 
total systematic errors on the SZE-derived 7^ value in A2204 
is roughly 2:1. 

The low systematic uncertainties in our analysis make it pos- 
sible to lower the error budget on the temperature profile signif- 
icantly by stacking the SZE signal of several relaxed clusters 
(Basu et al., in preparation). It is true that a very long exposure 
will drive down the systematic uncertainties associated with the 
X-ray background correction, and the remaining systematic un- 
certainties in the flux calibration are small (< 5% for Chandra, 
Vikhlinin et al. 2005). A precise comparison between the two 
gas temperature profiles, derived from joint SZE/X-ray analysis 
and X-ray spectroscopy, will be the most promising way to ob- 
servationally constrain gas clumping and non-LTE effects near a 
cluster's virial radius. 

We mention here the recent advances made by the X-ray 
spectral analysis method to constrain gas temperatures out to r2oo 
using the Suzaku experiment (Reiprich et al. 2009, George et al. 
2009, Bautz et al. 2009). This is due to the low level of particle 
background in the Suzaku orbit as compared to XMM-Newton 
and Chandra. For Abell 2204 Reiprich et al. (2009) have con- 
strained the gas temperature near r2oo at 4.49^', ^ keV, including 
both systematic and statistical errors. This is far superior to the 
current uncertainties in the APEX-SZ measurement. However, 
the extended PSF of Suzaku limits its ability to spectroscopi- 
cally measure the gas temperature out to the cluster virial radius 
to only low redshift (z < 0.2) massive clusters (most of which 
are too extended for single-frequency APEX-SZ measurement). 
This also makes modeling of the gas temperature at the inner 
radial bins difficult. Joint SZE/X-ray temperature modeling with 
XMM-Newton, Chandra or ROSAT data for X-ray surface bright- 
ness is therefore promising for the majority of clusters out to 
high redshifts. 

4.5. Direct comparison of de-projected pressure profile with 
parametric models 

Applying Abel's inversion technique to the SZE map produces 
an unbiased and non-parametric estimate of the cluster pressure 
profile for a spherically symmetric system. This can be used to 
compare the usability of different parametric models, needed to 
extract cluster properties like M200, from SZE or X-ray mea- 
surements made within ^500- For example, parametric extension 
is unavoidable while using interferometric measurements of the 
SZE signal in low and intermediate redshift clusters. 

The de-projected pressure profile assuming spherical sym- 
metry for A2204 is presented in Fig[8] error bars show the lcr 
statistical uncertainties in the SZE measurement. We have plot- 
ted the best-fit spherical isothermal /3-model and Nagai model 
(Nagai et al. 2007) fits on this profile. The fits are limited only 
to data within 6' radius, to mimic an SZE observation with lim- 
ited spatial extent that uses parametric model fitting to extrap- 
olate out to the cluster virial radius. We used the Nagai pro- 
file parameters as used by the recent SZA analysis of pressure 
profiles (Mroczkowski et al. 2009), with parameters (a, b, c) = 
(0.9,5.0,0.4). Our fitted scale radius is much larger than the 
predicted value of r p w r5oo/1.3 (7.3' in our fit), although it is 
strongly degenerate with the normalization factor. Similarly, the 
/3-model parameters [r c , f3] are also highly degenerate. We use fit 
values r c — 1 .4' and /3 = 0.5 1 for this comparison; setting /3 = 1 
we obtain r c = 2.0' which provides a marginally better fit to the 
peak SZE decrement in the deconvolved map. 

The yS-model is found to provide a poor extrapolated fit 
to the pressure profile even at r$oo (FigH), whereas the Nagai 
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Fig. 8. Comparison of the de-projected pressure profile with 
commonly used parametric models. The data points with errors 
show the result of direct inversion of the SZE temperature decre- 
ment map using Abel's integral. The red (solid) line is the best fit 
isothermal y6-model to the APEX-SZ map, fitted within 6' radius. 
The blue (dashed) line is the best fit Nagai model fitted within 
the same radius. The limited fitting radius was used to illustrate 
how well parametric modeling fitted within cluster central re- 
gions does reproduce the gas pressure in the outskirts. 



model provides a much better fit. The two outer bins represent- 
ing roughly the values at rgno and r2oo in our analysis have 3<x 
upper limits at 2.25 x 10 4 keV cirT 3 and 1 .01 X 10 4 keV cm~ 3 , re- 
spectively. This puts the values predicted by the Nagai model at 
these radii at roughly 1 .7cr above our measured values, whereas 
the y6-model predictions are at 4-5<x off. Note the limited spatial 
range used in this comparison study; a more accurate y6-model fit 
for the full SZE-derived pressure profile is possible given the de- 
generacy of the fitting parameters. 

We have also attempted to fit the de-projected pressure pro- 
file using a polytropic model for the gas, where the pressure and 
density are related by the relation P{r) cc n e {r) y . Many authors 
continue to use this model, e.g. Afshordi et al. (2007) use it to 
constrain cluster pressure profiles from WMAP data. This model 
is found to be too steep for the pressure profile near the cluster 
center, although in the outer regions (r > 1 Mpc) the fit is good 
(using y = 1.2, e.g. Finoguenovet al. 2001). A combination of /}- 
model in the cluster center and polytropic model in the outskirts 
can be used to fit the entire pressure profile, in particular to avoid 
the cuspiness of the Nagai model at the center. A more compre- 
hensive analysis of the different parametric models to describe 
the SZE-derived pressure profile near r2oo will be discussed in a 
future paper. 



5. Mass and entropy profiles of the ICM 

If the ICM is in hydrostatic equilibrium (HSE) within the DM 
gravitational potential, the gas temperature reflects directly the 
depth of the potential well. The ratio between gas mass and total 
mass as function of radius shows the amount of baryons that is 
contained in the ICM. A low value of the ICM mass fraction, or 
a falling gas entropy profile, can indicate the existence of multi- 
phase ICM with non-thermal pressure support near the cluster 
virial radius, and physical processes hitherto unexplored in nu- 
merical cluster simulations. 
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Fig. 9. Cumulative gas mass and total mass profiles in Abell 
2204 from the SZE/X-ray joint analysis. The black data points 
(boxes) show the result using the de-projected temperature and 
density profiles, assuming hydrostatic equilibrium condition. 
The green data points (triangles) show the results from XMM- 
Newton analysis by ZF08 under the same assumption. The blue 
solid line is the best fit NFW model from weak lensing analysis 
by Corless et al. (2009), the hatched region at r2oo indicates their 
quoted errors in M2oo- The lower red dashed line shows the gas 
mass profile computed directly from the de-projected electron 
density. 

5. 1 . Gas mass and total mass distribution 

The total mass, M tota i, is obtained by solving the hydrostatic 
equilibrium equation assuming spherical symmetry (e.g. Sarazin 
1988): 



M tota i(< r) = - 



k B T e (r) r 
Gymip 



dlnn e (r) d\nT e (r) 
dlnr dlnr 



(5) 



where T e (r) and n e (r) are the electron gas temperature and 
density radial profiles, and fj. = 0.62 is the mean molecu- 
lar weight per hydrogen atom, assuming primordial abundance. 
As can be seen, the total mass is primarily a function of the 
gas temperature, and only weakly dependent (through logarith- 
mic derivatives) on the slopes of the density and temperature 
profiles. Therefore, the uncertainties in our temperature mea- 
surement are directly reflected in the total mass profile. The 
gas mass is computed directly from the de-projected density 
as Pgas(r) = /J e m p n e (r), with pi e = 1.17 the mean molecular 
weight per electron. The gas mass fraction is simply the ratio: 
/ gas (< r) = M gas (< r)/M total (< r). 

Results from the non-parametric mass modeling are shown 
in Fig(9j in comparison with results from recent X-ray and weak 
lensing analyses. The mass profile is in excellent agreement with 
the X-ray results obtained under the same assumptions of spher- 
ical symmetry and hydrostatic equilibrium (ZF08). The solid 
line in Fig (prefers to the best fit NFW model for A2204 from 
weak lensing analysis by Corless et al. (2009). The mean mass 
profile indicates a slowly rising integrated mass near 7-500, re- 
sulting from the rapid fall in gas temperature near this radius, 
although the results are still consistent with the NFW model 
from Corless et al. (2009) within lcr errors. ZF08 reports a 
value of M 50 o = (5.8 ± 1.6) x 10 14 M Q from the X-ray model- 
ing, and the weak lensing analysis under a spherical prior gives 
M500 = 5.3 x 10 14 M e . Our SZE/X-ray joint de-projection anal- 
ysis predicts M500 = (3.7 + 1.2) x 10 14 M Q , somewhat lower than 
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Fig. 10. Ratio of the total cluster mass derived from the hydro- 
static equilibrium assumption, to the best fit NFW model ob- 
tained from the optical data. The cumulative total mass derived 
from APEX-SZ temperature measurements are consistent with 
the NFW model by Corless et al. (2009) under a spherical prior. 
For comparison, three data points from the stacking analysis for 
X-ray (XMM-Newton) to weak lensing (Subaru) mass ratio in 5 
low redshift relaxed clusters are shown (Zhang et al. 2010). 



the X-ray and weak lensing results. The lcr upper limit on M200 
from our analysis is 1.01 x 10 15 M , again consistent with the 
Corless et al. (2009) maximum likelihood model prediction of 
M 2 oo = 0.71+q^ x 10 15 M Q (shown by the blue hatched region 
in FigO. 

To show more clearly the deviation of the non-parametric 
mass modeling under the HSE assumption from the weak lens- 
ing mass, we plot in FigJTU] the ratio of the hydrostatic mass 
and the mass derived from the weak lensing analysis. Optical 
observations of A2204 are complicated by its low galactic lati- 
tude and presence of a My =5.6 star 4.3' away from the cen- 
ter, making its shear profile in the cluster center extremely noisy 
(Clowe & Schneider 2002). Thus we have used the profile from 
the NFW model fit to the weak lensing data instead. The ra- 
tio obtained is mostly consistent with 1 within lcr statistical er- 
rors. The HSE assumption is expected to under-estimate the total 
mass by 15-20% near the virial radius, due to the stochastic gas 
motions caused by infalling matter (Nagai et al. 2007, Ameglio 
et al. 2009). However, the current uncertainties on the APEX- 
SZ measurement of a single cluster are too large to confirm any 
such trend. For comparison, we have shown in Fig[l0]the recent 
results from the stacking analysis for the Mhse/Mvl ratio using 
XMM-Newton and Subaru data of 5 relaxed clusters (Zhang et al. 
2010). Note however, that this joint X-ray/weak-lensing stack- 
ing analysis uses the actual weak lensing shear measurements in 
clusters and not the best-fit NFW profiles, which should cause 
their ratio to be closer to unity than ours. 

The integrated baryon fraction of the ICM as function of ra- 
dius is computed directly by dividing the gas mass by the total 
mass obtained from the HSE assumption. The results from our 
non-parametric analysis are shown in FigQj] The cosmic baryon 
fraction obtained from the WMAP 5-year result (Q.i,/Q. m = 
0.165 ± 0.009, Dunkley et al. 2009) is shown in horizontal dot- 
dashed line. In the inner region of the cluster (r < rgoo) the gas- 
to-mass ratio is clearly much lower than the cosmic baryon frac- 
tion, and there is an indicative trend of increasing ICM mass 
fraction at larger radii. Near r2oo the cumulative value of / gas 
is statistically consistent with the cosmic value. Low values of 
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Fig. 11. The gas mass fraction obtained from the SZE-X/ray joint 
non-parametric analysis (red squares). The hatched region is the 
measurement from Afshordi et al. (2007) from the stacking anal- 
ysis of 193 massive clusters in the WMAP 3-year data. The solid 
line is the prediction from the isothermal yS-model fit to the SZE 
deconvolved map (r c = 1.78', (3 = 0.812), and the horizontal 
dot-dashed line represents the cosmic baryon fraction from the 
WMAP 5-year result (Dunkley et al. 2009). 



Fig. 12. Comparison of entropy profile in A2204 from the SZE- 
X-ray joint analysis, and Chandra X-ray analysis by Sanders et 
al. (2008). The red triangles are from the Chandra data, shown 
without errors. The dotted line is the best fit entropy profile of 
the form K(r) cc r • using the data points within 8' radius. The 
present statistical errors fail to give any definitive indication for 
low entropy gas near r2oo- 



gas mass fraction near cluster centers is well known from X-ray 
studies; Vikhlinin et al. (2006) have shown the value of / gas at 
7"2500 for a sample of nearby relaxed cluster to be in the range 
0.04 — 0.1. It is interesting to note, however, that the gas mass 
fraction in Abell 2204 at 7-2500 is significantly lower than previ- 
ous APEX-SZ measurements of this ratio in non-relaxed clus- 
ters. The integrated / gas values in Abell 2163 at r > ^500 were 
found to be consistent with the cosmic baryon fraction (NBP09), 
and HL09 measured the integrated / gas for the Bullet cluster (IE 
0657-56) within r 2 5oo and 1 .42 Mpc to be in the range 0. 1 8-0.22. 
One obvious explanation for the higher / gas value at the center of 
dynamically complex systems can be due to the fact that merg- 
ing activity will most likely cause the gas to remain at the cen- 
ter while the dark matter halos are separated, thus causing an 
increase in the gas-to-mass ratio. Also the central AGN in the 
strong cool core cluster A2204 can be responsible for driving 
out the gas from the innermost region (Bhattacharya et al. 2008, 
Puchwein et al. 2008). 

The / gas profile predicted from the isothermal /3-model (fit 
to the SZE data) is shown in FigJTT] As noted in NBP09, our 
non-parametric modeling shows clear departure from the typ- 
ical isothermal yS-model prediction of / gas — > at the cluster 
center. We have also shown for comparison the results for the 
stacking analysis of 193 massive clusters with Tx > 3 keV 
from the WMAP 3-year data by Afshordi et al. (2007). Note 
that the resolution of WMAP does not allow a direct measure- 
ment of the mean pressure profile in clusters down to 0.1r2oo, the 
hatched region is the prediction from their numerical simulations 
with P gas > prior that is consistent with the WMAP stacked 
measurement. The stacking signal from WMAP data predicts a 
higher ICM mass fraction than seen in A2204 near the center, 
although beyond r25oo they are consistent with each other within 
lcr. The difference near the center is expected, since the sample 
of Afshordi et al. (2007) contains both relaxed and non-relaxed 
clusters, which results in a higher average / gas value. If simi- 
lar low gas-to-mass ratios are found consistently from SZE ob- 
servations of massive relaxed clusters, then in parallel with the 
currently favored lower value of cr 8 parameter (Komatsu et al. 



2009), this will cause significantly low cluster yields in blind 
SZE surveys. 

5.2. The entropy profile of the ICM 

The entropy profile can be considered a more fundamental prop- 
erty for analyzing the thermodynamic state of the ICM than den- 
sity or temperature, as it depends directly on the history of heat- 
ing and cooling within the cluster. Consequently, cluster entropy 
profiles have received significant attention in X-ray studies (e.g. 
Voit et al. 2005, Morandi & Ettori 2007, Cavagnolo et al. 2009), 
but no direct measurement of entropy from SZE-derived tem- 
peratures have been done. The latter can potentially overcome 
the biases inherent in the X-ray spectral analysis, caused by gas 
clumping and substructures, and also from multiple temperature 
components near the cluster core. In this final part of our work 
we present the first SZE-derived entropy measurement in a clus- 
ter. 

We adopt the standard definition of gas entropy used in the 
X-ray literature: K - T e n e 1 (Ponman et al. 1999). This re- 
lates to the classical thermodynamic entropy in an ideal gas 
as s — lnA^ 2 + constant. Simulations for self similar cluster 
models predict an entropy profile in the form of a power law: 
K{r) oc r u (Tozzi & Norman 2001, Voit et al. 2005), except 
at the very core of the cluster where excess entropy due to non- 
gravitational heating processes tend to flatten the entropy profile. 
This general behavior has been verified from numerous X-ray 
observations (e.g. Lloyd-Davies et al. 2000, Morandi & Ettori 
2007, Cavagnolo et al. 2009). 

The entropy profile of Abell 2204 obtained from the SZE 
temperature measurements is shown in Fig [12] Also shown are 
the radial entropy values obtained by Sanders et al. (2009) from 
the Chandra measurements (triangles). The statistical errors in 
their measurements are comparable with the symbol sizes and 
much smaller than the present SZE measurement errors (al- 
though note that the systematic uncertainties should dominate 
the errors in the outer Chandra bins, which was not shown by 
Sanders et al. 2009). The two measurements agree within the lcr 
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results between two clusters. A fit with the expected power law 
scaling excluding the inner- and outermost bins for A2204 gives 
* 2 /d.o.f. = 7.4/5, and for A2163 it is x 2 /d.o.f. = 3.2/5. 

To conclude, we are able to demonstrate the correlation of 
the "entropy floor" with the dynamical state of a cluster (FigQjJ) 
for the first time using SZE imaging data, independent of X-ray 
spectral analysis. This correspondence shows the potential for 
the current de-projection technique, using X-ray surface bright- 
ness maps and already available SZE imaging data from multi- 
pixel bolometer array experiments, to select relaxed clusters 
from a large sample where deep X-ray observation is not avail- 
able/required. This will be important, for example, for selecting 
relaxed clusters to constrain cosmology with the ICM baryon 
fraction. 



Fig. 13. Comparison between entropy profiles for A2204 (black, 
diamonds) and A2163 (blue, squares). The increased entropy 
value in the central region of A2163 supports the merging na- 
ture of this cluster, whereas the continually decreasing entropy 
towards the center in A2204 shows its dynamically relaxed state. 

uncertainties of our analysis. The agreement between the APEX- 
SZ and Chandra values within the central 1' radius shows that 
any downward bias on the APEX-SZ value due to PSF smearing 
is sufficiently small. The dotted line in Fig [12] is the power-law 
prediction from self-similar cluster models, fitted to our entropy 
measurement within 8' radius. 

A flattening of the entropy profile near the virial radius of 
clusters has been measured only very recently from X-ray spec- 
tral analysis with Suzaku data: for PKS 0745-191 (z = 0.10, 
George et al. 2009) and Abell 1795 (z = 0.06, Bautz et al. 2009). 
Note that the result of Bautz et al. (2009) was obtained by ex- 
trapolating the values within rsoo- Such flattening or dropping 
entropy profile is another indicator of the non-thermal pressure 
support in the cluster outskirts. From the mean de-projected ra- 
dial density and temperature profiles, we also see similar flatten- 
ing for Abell 2204 near rsoo (FigfT2l. but the statistical uncer- 
tainties in these SZE measurements are too large which makes 
our profile compatible with the power law scaling of the self 
similar cluster models. 

The 1 arcmin resolution of the APEX-SZ experiment is suffi- 
cient to compare the gas entropy at the cluster cores for low and 
intermediate redshift massive clusters (r2oo ^ 8')- As mentioned 
before, this comparison is important in analyzing the cluster 
dynamical states, by comparing the extent of non-gravitational 
heating or cooling at the cluster core. We have used the results 
for Abell 2163 from the analysis in NBP09, which is a good ex- 
ample of a dynamically complex cluster which most likely has 
undergone a merging event (Maurogordato et al. 2008). The data 
for A2163 is analyzed with the same radial binning as for A2204 
used in the current paper. The comparison of the entropy profiles 
between these two clusters is shown in FigQj] The difference 
in the central entropy values is clearly seen. Although the cen- 
tral bin errors are possibly under-estimated by a small amount 
because of neglecting PSF extension, the statistical significance 
of the entropy difference (~ 6cr) is sufficiently high. Outside the 
central ~ 2 arcmin the two profiles are consistent with each other 
and both follow the oc r 11 scaling law. It is interesting to note 
that the merging system A2163 shows better agreement with the 
power law scaling than the relaxed cluster A2204. Using a x 2 
statistic to describe the goodness of fit is problematic because of 
the noise correlation, nevertheless, it can be used to compare the 



6. Conclusions 

1. We describe the detailed application of a direct, nonpara- 
metric de-projection method of cluster density and tempera- 
ture profiles, using APEX-SZ and XMM-Newton data. The 
method was presented in NBP09, the current paper builds 
upon the previous work by applying this technique to the 
well-studied relaxed cluster Abell 2204. 

2. Analysis of both SZE and X-ray data have been done from 
their raw data sets, to create images and radial profile. In 
particular, we describe the creation of a set of SZE radial 
profiles, all consistent with the APEX-SZ measurement, to 
characterize the statistical uncertainties on the bin values 
and minimize the numerical errors in Abel's de-projection 
method. Our final results are dominated by the statistical un- 
certainties in the SZE data, the signal at r^oo is essentially an 
upper limit for A2204. We describe the different sources of 
systematic uncertainties and include them in the analysis. 

3. The decreasing gas temperature in the cluster outskirts is 
demonstrated for the first time from SZE measurements, us- 
ing a broad re-binning of the APEX-SZ (and X-ray) data. 
The temperature drop can be confirmed to 98% confidence 
level. We also compare the direct de-projected pressure pro- 
file with some parametric models, and show that the Nagai 
profile is adequate for modeling the gas pressure, within the 
current statistical uncertainties in APEX-SZ imaging of a 
single cluster. 

4. We re-perform the X-ray spectral analysis for the ICM tem- 
perature profile from publicly available Chandra data, pri- 
marily to find the changes from the recent calibration up- 
dates (CALDB 4), but also to show the effect of systematic 
uncertainties due to the background modeling in the X-ray 
spectral analysis. A comparison with the projected temper- 
ature profile obtained from SZE data confirms that our SZE 
derived temperature values are much less affected by system- 
atic uncertainties at large radii, in comparison with Chandra 
and XMM-Newton. Precise comparison between the SZE and 
X-ray spectroscopic measurements of the gas temperature in 
the cluster outskirts will be a promising method to constrain 
gas clumping and non-LTE effects. 

5. The integrated total mass profile is computed assuming hy- 
drostatic equilibrium for the cluster gas. The mass profile is 
in excellent agreement with the recent X-ray and weak lens- 
ing analyses. Our model prediction for M500 is (2.6 + 2.2) x 
10 14 hr l M Q . This is somewhat lower than the X-ray and lens- 
ing results but consistent within lcr errors. The upper limit 
on M200 from our analysis is in good agreement with the pub- 
lished NFW model fit from weak lensing analysis of A2204. 
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6. The ICM mass fraction as function of radius is computed us- 
ing the non-parametric modeling, and found to be below 0.1 
within 7-500. The low / gas value in A2204 in the cluster center 
is in contrast with the previous APEX-SZ measurement of 
this ratio for Abell 2163 and the Bullet cluster. 

7. We compute the ICM entropy profile from SZE/X-ray joint 
analysis and confirm the general agreement with the self- 
similar cluster model predictions within the present statisti- 
cal uncertainties. The significance of the APEX-SZ measure- 
ment of A2204 is not sufficiently high at r > rsoa to constrain 
the slope of the entropy profile in the cluster outskirts. 

8. We compare the entropy profiles of Abell 2204 (relaxed) 
and Abell 2163 (merging system), using the same non- 
parametric SZE/X-ray de-projection and radial binning, and 
find a clear entropy difference in their central 200 kpc. This 
corresponds to the different dynamical states of these two 
clusters and seen for the first time from SZE derived T gas 
measurement. 
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